Energetics of positron states trapped at vacancies in solids 
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We report a computational first-principles study of positron trapping at vacancy defects in metals 
and semiconductors. The main emphasis is on the energetics of the trapping process including the 
interplay between the positron state and the defect's ionic structure and on the ensuing annihilation 
characteristics of the trapped state. For vacancies in covalent semiconductors the ion relaxation 
is a crucial part of the positron trapping process enabling the localization of the positron state. 
However, positron trapping does not strongly affect the characteristic features of the electronic 
structure, e.g., the ionization levels change only moderately. Also in the case of metal vacancies 
the positron-induced ion relaxation has a noticeable effect on the calculated positron lifetime and 
momentum distribution of annihilating electron-positron pairs. 
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I. INTRODUCTION 

Positron annihilation spectroscopy^ is widely used in 
studying vacancy-type defects in metals and semiconduc- 
tors. This is based on the fact that vacancies, in neutral 
or negative charge states, act as efficient positron traps 
due to the reduced repulsion of positive ions. Also the 
electron density is reduced at vacancies and by measur- 
ing the ensuing lifetime increase with respect to the de- 
localized positron bulk state one can estimate the extent 
of the open volume in defects. The momentum distri- 
bution of annihilating electron-positron pairs measured 
by Doppler broadening spectroscopy is at high momenta 
specific for the annihilating core electrons and thereby it 
enables the chemical identification of atoms surrounding 
vacancies. 

For a given sample, the identification of the most 
abundant open- volume defect type, such as vacancy 
and vacancy agglomerates or vacancy-impurity complex, 
is based on the knowledge of the general behavior of 
positron annihilation results, i.e., on the information how 
the measured annihilation characteristics depend on the 
open volume or the chemical environment of defects. 
This knowledge has been acquired by measuring well- 
characterized reference samples including well-annealed 
(perfect bulk) materials as well as defected materials such 
as those containing monovacancies due to electron irradi- 
ation. Moreover, theoretical predictions of positron an- 
nihilation characteristics significantly support the defect 
identification.^ 

In order to interpret experimental results it is also im- 
portant to understand the positron trapping process in 
detail. The trapping coefficient is an important quan- 
tity in determining defect concentrations. Its values have 
been estimated also theoretically for model systems by 
assuming that the positron gives in the trapping pro- 
cess its binding energy to excited electron-hole pairs or 
to phonons."^'^'^ In the present computational work our 
main theme is the effect of the positron on the ion lattice 
during the trapping process. Because the ion vibration 
frequency, which is of the order of the Debye phonon fre- 
quency, is much larger than the positron annihilation rate 
the ions around the defect have time prior to positron 



annihilation to relax to minimize the total energy of the 
defect-positron system. The ion relaxation affects the lo- 
calization of the positron state and the annihilation char- 
acteristics. Below we argue using first-principles total- 
energy calculations that the positron-induced lattice re- 
laxation is indispensable for the existence of localized 
positron states at vacancies in covalent semiconductors 
and thereby it completely determines the ensuing an- 
nihilation characteristics. In these systems, the strong 
infiucnce of the positron is possible because the energy 
landscape as a function of the ion positions around the 
vacancy is very flat. Actually, the effect of the trapped 
positron is found so strong that it practically cancels the 
possible symmetry-lowering Jahn- Teller distortion of the 
vacancy. For vacancies in metals the influence of the 
trapped positron on the ion positions and especially on 
the energetics of the trapping process is smaller. How- 
ever, the trapped positron causes a small increase in the 
open vacancy volume and thereby noticeable chances in 
the positron lifetime and in the momentum distribution 
of the annihilating pairs. It is important to note this 
from the modeling point of view. 

In the case of semiconductors, positron annihila- 
tion has been used also to extract detailed informa- 
tion about the electronic structures of the defects, i.e., 
to determine the ionization levels of vacancy-type de- 
fects^iiiSiiSiiiS'iiii^'i^ or just to probe the changes in their 
charge statesJ^ In these experiments the charge state 
of the defect changes due to the thermal ionization or 
by illumination with light. The charge state change is 
then seen as a change in the positron lifetime reflecting 
electronic-structure-induced changes in the ion positions 
or as a dramatic change in the positron trapping rate 
when positive defects do not trap positrons. The ques- 
tion rising immediately in the first case is how much the 
positron-induced ionic relaxation affects the positions of 
the ionization levels by modifying or eventually breaking 
bonds between ions next to the vacancy. Our prediction 
is that although the changes in the ionic structure are 
rather large their effect on the ionization levels is minor. 

Finally, we would like to point out that the com- 
parison of calculated positron annihilation characteris- 
tics with the measured ones constitutes the fundamen- 
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tal test for theories describing electronic properties of 
materials and the positron-electron interactions as well 
as for many computational approximations. For delo- 
calized positron states in perfect bulk solids there exist 
several systematic comparisonai^ii^ii^ii^ but for positron 
states trapped at vacancy defects comparisons treating 
several materials and systems on the same footing are 
scarce. The reason may be in difficulties arising in the 
theoretical description, e.g. in the density-functional the- 
Qj, yi9,20 (^Qp^p) local-density approximation (LDA) 
for the electron exchange and correlation underestimates 
the energy band gap in semiconductors which may have 
severe consequences on the localized electron states and 
the ionic structure at defects. Moreover, the electron- 
positron correlation effects are known worse for localized 
positron states than for delocalized ones. Lastly, the bro- 
ken translational symmetry leads to computational ap- 
proximations such as the supercell method which requires 
large computer resources in order to show convergence of 
the results with respect to the supercell size. The aim of 
the present study is to remedy the situation by providing 
results for a representative set of materials. We consider 
metals with different lattice structures (close-packed Al, 
Cu, Mg, and body-centered-cubic Fe) and elemental (Si, 
Ge) and compound semiconductors (GaAs, GaN) with 
different degrees of bond ionicity. The structure of the 
present paper is as follows. In Sec. |TT] we review shortly 
the theory and computational methods used. The results 
are given and described in Sec. lIIII and they are discussed 
along comparisons with experimental results in Sec. IIVI 
Section IVl presents our conclusions. 



equation for the positron state ip+{r) reads as 

^+(r) = 0(r) + yeorr(ri-(r)), (1) 

where (/)(r) is the Coulomb potential due to electrons and 
nuclei, n_(r) the electron density and T4orr('^-) is the 
electron-positron correlation energy^! for a positron in a 
homogeneous electron gas with density n_. It has been 
shown that the effects of the above two approximations 
largely cancel each other's effects so that the conven- 
tional scheme and 2CDFT results for positron annihi- 
lation characteristics, for the total energy of the defect- 
positron system and for the positron trapping energy are 
very similari^ii^^ Besides due to the conceptual simpli- 
fication, we prefer the conventional scheme also because 
the 2CDFT functionals for finite positron densities are 
not accurately known. 

When we relax the ions surrounding a defect with a 
trapped positron we minimize the total energy which in 
the conventional scheme is the sum of the DFT total 
energy for the electron-ion system and the positron en- 
ergy eigenvalue. Thus, although our calculation is not a 
self-consistent 2CDFT calculation the positron state and 
the electron density are coupled via the ionic structure. 
In practice, we calculate the positron-induced Hellman- 
Feynman forces on ions using the so-called atomic super- 
position method (for details, see Ref. 18). 

The total annihilation rate A which is the inverse of the 
positron lifetime t is obtained from the overlap integral of 
the electron density n_(r) and positron density n-^.{r) = 
lV'+(r)P 



II. THEORY AND COMPUTATIONAL DETAILS 
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A. Theoretical models 

We perform first-principles electronic-structure calcu- 
lations based on DFT for various vacancy defects in met- 
als and in semiconductors. These calculations give the 
ionic positions by requiring that the total energy is min- 
imized. This is equivalent to vanishing of the Hellman- 
Feynman forces on ions, calculated from the ground- 
state electron density. The trapped positron state at 
a defect can be included by generalizing to the two- 
component density- functional theory^^ (2CDFT). For de- 
fects in semiconductors calculations optimizing the elec- 
tronic and ionic structures as well as the positron density 
within the 2CDFT have appeared,^i^i2ii25,26 

In the present work we apply the so-called conventional 
scheme in which (i) the localized positron density does 
not directly affect the average DFT electron density (the 
positron and its screening electron cloud form a neutral 
quasiparticle entering the system) and {ii) the positron 
state and annihilation characteristics are calculated in 
the LDA and at the zero-positron-density limit of the 
electron-positron correlation functionals. For example, 
this means that the potential entering the single-particle 



Above, 7(r7-_) is the enhancement factor for a positron a 
homogeneous electron gas^i with density n_ , and and 
c are the classical electron radius and the speed of light, 
respectively. We calculate the momentum distribution 
p(jp) of the annihilating electron-positron pairs using the 
state-dependent enhancement schemej^li^ I.e., 



(3) 



where ipj (r) and 7j are the wave function and the state- 
dependent enhancement factor (in the LDA) for the elec- 
tron state j. 

Besides the agreement with the 2CDFT the results ob- 
tained with the conventional scheme compare reasonably 
well also with experiment.— First and foremost, the 
measured changes in the positron lifetime and in the rel- 
ative changes in the core annihilation rate between the 
vacancy and bulk states are reproduced. 

The most important aspect of the present work is the 
energetics of the defect-positron system. We define the 
positron trapping energy at a vacancy defect (the en- 
ergy released in the trapping process) as the total energy 
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difference between the systems of (z) a defect and a de- 
localized positron and (m) the same defect trapping a 
positron. Within the conventional scheme we obtain 

Et = Ai^tot ^iE + e+ - (E,, + £+ 

= (4uik-4ofcct)-(^e+-i?), (4) 

where E is the total energy of the electron-ion system of 
the defect supercell without the localized positron, £'e+ 
that with the localized positron, e^^^^. and e^cfcct 
energy eigenvalues of the positron in the delocalized bulk 
state and in the localized state at the vacancy, respec- 
tively. The last form in Eq. ^ shows that the trapping 
energy consists of the decrease of the positron energy 
eigenvalue and the increase in the (strain) energy stored 
in the ion lattice. In general, the Kohn-Sham eigenval- 
ues in the 2CDFT have no physical meaning but as we 
have only one positron in the lattice and we use the con- 
ventional scheme the electron-positron interactions affect 
only the positron energy eigenvalue e+ and the above 
analysis is justified. 

The ionization level e{Q/Q') between the charge states 
Q and Q' of a defect is defined as the position of the 
chemical potential fj,e relative to the top of the valence 
band Ey so that the total energies of these two charge 
states are equal. I.e., we solve for in 

E2,t + Q{Ey + fie) = i?tot + Q'iEv + Me), (5) 

where i?^^. is the total energy of the supercell with the 
defect in the charge state Q and the term Q{Ey + fie) 
arises because Q electrons are added to (Q > 0) or taken 
from (Q < 0) the electron reservoir at the chemical po- 
tential level Ey + fie- 

B. Calculation methods 

In the supercell approach we use one has to take care 
of well-known artifacts. First, the energy eigenvalue 
^buik delocalized positron in Eq. Q and the va- 

lence band maximum Ey in Eq. ([5]) are taken from the 
calculation for the perfect periodic bulk material. Be- 
cause the energy zeroes differ between different super- 
cells we align the effective potentials for the defect and 
bulk supercells far from the defect both in the case of 
electrons and the positron. Second, in order to avoid 
long-range Coulomb interactions between charged super- 
cells in the superlattice we use a neutralizing uniform 
background charge. The unphysical energy terms due 
to the monopole-monopole interactions between the pe- 
riodic images of the defects are corrected by the method 
by Makov and Payne. However, in the case of defects in 
GaN we do not apply these corrections because they lead 
to an overcorrection. Since we are primarily interested 
not in the absolute values of the ionization levels but in 
their changes due to the localization of the positron the 
energy corrections are actually not of utmost importance. 



Our computational methods are described in more de- 
tail in Ref. |18i and here we will give only the main 
features. We perform electronic structure calcula- 
tions within the LDA (Ref. Is^ ). The description of 
the electron-ion interaction is based on the projector 
augmented-wave (PAW) method^ implemented in the 
plane- wave code VASP (Refs. HHHIll). Using the PAW 
total charge density including the free atom core elec- 
trons the positron potential is constructed according to 
Eq. ^ and the lowest-energy positron state is calculated 
on a three-dimensional real-space point grid. 

The lattice constants of the perfect bulk lattices are 
optimized and used in the defect calculations to define 
the supercell volume. For Si, Ge, and GaAs we use cubic 
216-atom supercells. The Brillouin zone is sampled in the 
case of Si and GaAs using a 2'^ Monkhorst-Pack (MP) k 
point meshes^ whereas for Ge we use the L-point sam- 
pling in order to avoid the artificial hybridization of the 
deep level and band states in the LDA (see Ref. ^) . Cut- 
off energies for Si, Ge, and GaAs are 246 eV, 270 eV, and 
209 eV, respectively. We model wurtzite GaN using an 
orthorhombic 96-atom supercell, a 3^ MP k-mesh and 
a cutoff energy of 400 eV. The hexagonal close-packed 
(hep) Mg is modeled using an orthorhombic 48-atom su- 
percell, a 8 X 6 X 6 k-mesh and the cutoff of 263 eV. For 
the body-centered cubic (bee) Fe we calculate the mag- 
netic ground states using a cubic 54-atom supercell, a 8'^ 
MP k-mesh and a cutoff energy of 268 eV. For defects 
in the face-centered cubic (fee) metals Cu and Al we use 
cubic 108-atom supercells, a 6'^ and 8^ MP k-meshes and 
the cutoff energies of 342 eV and 301 eV, respectively. 



III. RESULTS 
A. Ion relaxation in positron trapping at vacancies 

We study first the energetics of the trapping process, 
i.e., the interplay between the lowering of the positron 
energy eigenvalue and the energy stored in the strained 
lattice around the vacancy. Vacancies in Al and in Si rep- 
resent two very different types of behaviors. In order to 
facilitate the scanning of the energy landscape when the 
ions relax due to the positron-induced forces we consider 
only the breathing-type ion relaxation. Thus, the point 
symmetries of the Al and Si vacancies are constrained 
to be Oh and T^, respectively. The reaction coordinate 
is the relaxation of the nearest-neighbor ions of the va- 
cancy from their ideal lattice positions. The positions of 
the other ions in the supercell are optimized. 

The results for the Al vacancy are shown in Fig. [ija). 
The energy of the electron-ion system (the uppermost 
curve) and the positron energy eigenvalue (the lowest 
curve) as well as their sum, the total energy of the sys- 
tem (the curve in the middle), are shown as a function of 
the relaxation of the nearest-neighbor ions. The energy 
zero is chosen to be the total energy of the vacancy and 
the trapped positron with ions relaxed without positron- 
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FIG. 1; Configuration-coordinate diagrams for (a) the mono- 
vacancy in Al and (b) the neutral monovacancy in Si. The 
dotted fines show the positron energy eigenvalue (relative to 
the one in perfect bulk), dash lines correspond to the energy 
of the lattice and the delocalized positron and solid lines the 
total energy of the defect-positron system as functions of the 
relaxation of the vacancy. Positive (negative) sign denotes 
outward (inward) relaxation. The points A, B, and C denote 
different stages in the positron trapping process. 



induced forces. Then the uppermost curve corresponds 
also to the total energy of the vacancy and a delocalized 
positron. The smallest relaxation shown corresponds to 
the equilibrium ion positions of the Al vacancy without 
the trapped positron. The picture of the positron trap- 
ping process is clear. First, a fast (fast compared to the 
time scale of ionic movement) vertical Franck-Condon 
shift of the positron from the delocalized bulk state to 
the localized ground state at the vacancy takes place via 
electron- hole excitation [A— >B in Fig. [Ha)]. Then, the 
ions move slightly outward to minimize the total energy 
of the vacancy-positron system (B— >C). The escape of the 
positron from the trapped state via thermal processes, 
the so-called detrapping process, is very unlike due to 
large separation of two uppermost total energy curves. 

The energetics for the neutral Si vacancy trapping a 
positron is depicted in Fig. [Ijb). Actually, in this case 
we cannot find a bound positron state for the strongest 
inward relaxations and therefore the curves join at the 
leftmost point corresponding to the symmetric equi- 
librium relaxation without the trapped positron. At 
smaller inward relaxations a bound state exists and the 
positron energy eigenvalue decreases. Surprisingly, the 
energy stored in the electron-ion system and the lower- 



ing of the positron energy eigenvalue cancel each other 
rather accurately and the sum curve is very flat over a 
large range of ionic relaxation. The total energy shows a 
maximum and at small outward relaxations a minimum 
which gives a trapping energy [Eq. ([4])] of 0.05 eV. In the 
trapping process a thermalized positron would "clear" a 
larger empty volume and a slightly deeper potential well 
for itself. In this way the trapping process is analogous 
to the self-trapping of electrons or holes in small polaron 
states in ionic crystals (see for example Ref. IstI ) . The flat 
total energy surface means that at finite temperatures the 
ions can be quite far from their absolute minimum energy 
configuration. However, the entropy contribution to the 
free energy would favor a larger open volume; The situ- 
ation that all the four nearest-neighbor atoms are very 
close each other (strong inward relaxation) is very un- 
likely due to the small corresponding phase-space volume. 
As a consequence (assuming that positron detrapping is 
a vertical Franck-Condon process) there would also be 
a finite effective detrapping energy related to the dis- 
tance between the two uppermost total energy curves in 
Fig. [ijb). The detrapping energy would be of the order 
of tenths of an eV. 

In conclusion, the main characteristic differences be- 
tween positron trapping at Al and Si vacancies are the 
much larger change in the positron energy eigenvalue for 
the Al vacancy and the much larger ion relaxation at the 
Si vacancy. The repulsive effect of the localized positron 
is stronger than one might expect on the basis of its small 
charge. The zero-point motion of the positron increases 
the force on the neighboring ions because the positron 
density penetrates closer to their nuclei. For example, 
we estimate that in the case of the Si vacancy the force 
due to a classical positive unit point charge at the cen- 
ter of the vacancy is only ~ 50% of the force due to the 
localized positron. For the Al vacancy the Oh symme- 
try persists also without constraints but for the neutral 
Si vacancy a symmetry lowering Jahn- Teller distortion is 
expected. In fact, when the defect is relaxed without a 
localized positron we find a Jahn- Teller distortion with 
the D2d symmetry that lowers our calculated trapping 
energy slightly so that it even becomes negative. When 
the positron is trapped at the vacancy its repulsion prac- 
tically restores the Td symmetry of the vacancy. Within 
the numerical accuracy we can consider the trapping en- 
ergy to be practically zero. The flatness of the energy 
landscape will be a general characteristic feature of the 
Si vacancy trapping a positron. 



B. Trapped positron states and annihilation at 
vacancies 

Next we present our flrst-principles results for various 
semiconductors and metals. Figure [2] shows the calcu- 
lated positron densities at Al, Fe, and Si vacancies and 
in corresponding defect-free lattices. Thus, examples of 
fee and bcc metals and tetrahedrally-bonded semicon- 
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FIG. 2: Positron densities in monovacancies (left, contour 
spacing 0.01 A~^) in Al, Fe and Si and corresponding perfect 
lattices (right, contour spacing one tenth of the maximum 
value). The dots in the figures denote the locations of the 
nuclei on the plane. 

ductors are considered. The vacancies in metals localize 
the positron state effectively whereas in Si the positron 
density tends to leak along the open interstitial channels, 
which is reflected also in the smaller maximum value of 
the positron density. Because of the higher coordination 
number the Coulomb repulsion due to the nuclei is larger 
in the interstitial regions in the fee and bcc lattices than 
in the open interstitial channels in the diamond struc- 
ture of Si. Therefore the positron energy eigenvalue will 
decrease in the trapping process more in the fee and bcc 
metals than in semiconductors 

More quantitatively, we calculate positron trapping en- 
ergies at vacancies and analyze how localized positrons 
affect the volumes and symmetries of the defects. More- 
over, for vacancies in semiconductors we determine ther- 
modynamical ionization levels in the energy band gap 
with and without a trapped positron. Our main results 
are presented in Table |T] and they are discussed in the 
following subsections. 

1. Trapping energies 

The positron trapping energies with components giv- 
ing the lowering of the positron energy eigenvalue and 
the energy stored in the ion lattice are shown on the 
leftmost data columns in Table |T] [see the last form in 
Eq. (ID)]. Several trends can be seen. First of all, the 



positron trapping energy at vacancies (V) in metals is 
typically clearly larger than that in semiconductors, i.e., 
the values for Vai, Vpc, and Vcu are of the order of 2 eV 
whereas the values for Vsi and Voa in GaAs are at most 
a few tenths of an eV. As can be seen in Table [T] this 
difference originates mainly from the fact that the lower- 
ing of the positron energy eigenvalue is larger in metals 
than in semiconductors. This in turn reflects the reduc- 
tion of the positron-nucleus Coulomb repulsion which is 
larger when a vacancy is created in metal lattices with 
a larger atomic density and higher coordination number 
than in tetrahedrally coordinated semiconductor lattices 
with large open interstitial channels. However, there are 
exceptions from this general trend. Among the vacancies 
in metals, VMg has a very low positron trapping energy 
which reflects the relatively low atom and electron densi- 
ties. Vca in GaN has a positron trapping energy similar 
to metals. This is no wonder because due to the size dif- 
ference between the Ca and N atoms the Ca atom density 
in GaN is more than 60% of that in the Ga metal. The 
high atom density increases the positron-nucleus repul- 
sion in the perfect GaN lattice and the lowering of the 
energy eigenvalue in trapping. 

According to Table [1] the energy stored in to the lattice 
relaxation during the positron trapping process is of the 
same order of magnitude, ~ 0.5 eV. However, the crucial 
difference between typical metals and semiconductors is 
that for metals this energy is only a fraction of the energy 
which the positron gains in the lowering of the energy 
eigenvalue while for semiconductors it is of the same order 
of magnitude. Indeed, in the case of Vsi, and Vca in GaAs 
the two components of the trapping energy nearly cancel 
each other. 

Our calculations predict that the trapped positron 
state and the accompanying ionic relaxation at Vsi is 
only a metastable configuration and that the global en- 
ergy minimum corresponds to an unperturbed vacancy 
and a delocalized positron. The energy barrier between 
these two minima is only some tenths of meV larger than 
the absolute value of the predicted negative trapping en- 
ergy. The actual situation differs from that depicted in 
Fig. [TJb) because allowing the symmetry-breaking the 
Jahn- Teller effect lowers the energy of the vacancy re- 
laxed without the trapped positron. 

In the case of Vqc and Vas in GaAs we find a bound 
positron state when the ions are frozen at ideal lattice 
positions or the nearest neighbor ions are relaxed out- 
wards. But when we start optimizing the ion positions 
the vacancy relaxes strongly inward destroying the bound 
state. Vgg can be contrasted with Vsi; The increase in 
the lattice constant from Si to Ge does not compensate 
the increase in the ion size and not even a metastable 
configuration with a bound positron state is found. For 
Vn in GaN we did not found a bound positron state even 
when the ions neighboring the vacancy were frozen at po- 
sitions corresponding to reasonable outward relaxations. 
In the latter case this can be explained by the small size 
of the N ion. 
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TABLE I: Positron trapping energies and their decompositions into the decrease of the positron energy eigenvalue and increase 
of the energy of the lattice, ionization levels, relative volume changes and resulting point symmetry groups, positron lifetimes r 
and relative W parameters for various vacancy defects in different charge states Q in bulk solids. Results are calculated with a 
positron trapped at the defect except for the ones in parenthesis which are obtained without the trapped positron. Ionization 
levels are given with respect to the valence band maximum. The relative changes {V — Vo)/Vo in the vacancy open volume are 
calculated from the volumes of the polyhedra defined by the nearest-neighbor atoms in the ideal (Vb) and relaxed {V) lattice 
positions. Negative (positive) values correspond to inward (outward) relaxation of the nearest-neighbor atoms. The positron 
lifetimes at defects should be contrasted with our corresponding computational lifetimes for the bulk solids which are Si: 208, 
Ge: 213, GaAs: 212, GaN: 131, Fe: 87, Mg: 219, Cu: 95, and Al: 159 ps. 



Defect 


Q 


Et 


e+ - =■+ 

= bulk ^defect 




e(Q/Q"l) 


iV~Vo)/Vo 


Symmetry 


T 








(eV) 




(eV ) 


(eV) 


(eV) 


(%) 




(ps) 




V in Si 





-0.17 




u.oy 


0.86 


0.47 (0.47) 


+20.5 (-43.1) 


Td (I>2d) 


260 


0.55 




-1 


-0.17 




n fin 
U.oU 


0.97 




(-54.8)" 


Ta (Dsd)" 


255 


0.57 


Vbe in Ge 







(no 


minimum 


with bound e+) 


(0.05) 


(-42.5) 


(^2d) 








-1 




(no 


minimum 


with bound e^) 




(-45.3) 


P2) 






Vba in GaAs 


-2 


0.39 




0.75 


0.36 


0.74 (0.81) 


-13.0 (-37.3) 




237 


0.68 




-3 


0.46 




0.84 


0.38 




-16.2 (-37.5) 




234 


0.69 


Vas in GaAs 







(no 


minimum 


with bound e"^) 


(0.21) 


(-41.1) 


(^2d) 








-1 




(no 


minimum 


with bound e^) 




(-51.1) 


{D2,) 






Vba in GaN 


-2 


1.82 




2.05 


0.23 


1.22 (1.43) 


-f 59.0 (4-29.7) 


Czv (Cst,) 


216 


0.57 




-3 


2.04 




2.28 


0.24 




4-60.7 (+29.4) 




216 


0.55 


Vn in GaN 







(no 


bound e^ 


state) 


(2.62) 


(-7.7) 


(Cs.) 








-1 




(no 


bound e+ 


state) 




(-19.5) 


{Czv) 






V in bcc Fe 




1.67 




2.34 


0.67 




+8.4 (-6.2) 


On (Oh) 


159 


0.75 


V in hep Mg 




0.34 




0.49 


0.15 




+6.3 (-2.9) 


Csh (Csh) 


289 


0.56 


V in fee Cu 




2.20 




2.58 


0.38 




+7.4 (-3.8) 


Oh (Oh) 


163 


0.74 


V in fee Al 




1.89 




2.35 


0.46 




+8.8 (-5.1) 


Oh (Oh) 


242 


0.78 



"Split vacancy. 



2. Defect geometries 

The relative changes {V — Vq)/Vq in the vacancy open 
volume are also given in Table HI Here Vb and V refer 
to the volumes of the ideal (atoms at the ideal lattice 
sites) and relaxed vacancies, respectively. Vq and V are 
calculated as the volumes of polyhedra restricted by the 
nearest-neighbor atoms of the vacancy. The numbers in 
the parenthesis show that without trapped positrons va- 
cancies in metals and in typical semiconductors have a 
tendency to shrink, i.e., the nearest neighbor atoms relax 
inward toward the center of the vacancy. The relaxation 
is large for semiconductors whereas in metals the ions re- 
main close to their ideal lattice positions. In GaN the N 
atoms neighboring Vca relax outward which reflects again 
the role of large Ga atoms in determining the lattice con- 
stant. The trapped positron increases the open volume. 
The effect is very dramatic in the case of semiconductors 
in which the volume increase of the vacancies is several 
tens of percents. For the neutral Vsi it is even of the 
order of 60%. The changes in the vacancy open volume 
are reflected in the positron lifetime and the momentum 
density of the annihilating electron-positron pairs (see 
Ref. [isl ). The effect is strong also in the case of metals 
although the changes in the ionic relaxations are smaller. 



The strong effect of the trapped positron is seen also 
in the results for the vacancy point symmetry in semi- 
conductors. Fig. [3]Ja) illustrates the case of the neutral 
Vsi- Without the trapped positron the Jahn- Teller ef- 
fect lowers the point symmetry to D2d and we see that 
the dangling bonds pointing toward the center of the 
vacancy hybridize to two pairs of bonds between the 
nearest-neighbor atoms. In Table [T] the symmetry of the 
plain singly negative Viji is D^d corresponding to the split- 
vacancy configuration where one of the atoms neighbor- 
ing the vacancy relaxes so that a divacancy with an atom 
in the center results. These results for the Si vacancy are 
in good agreement with previous LDA results j"^^'^^ With 
a positron localized at Vsi the ideal lattice point symme- 
try Td is practically restored. As seen in Fig.[3]Jb) for the 
neutral Vsi the strong positron repulsion increases the 
distances and weakens the bonds between the nearest- 
neighbor atoms of the vacancy and the atoms end up at 
the ideal-lattice point-symmetry positions within the nu- 
merical accuracy. However, the deep localized electron 
state is not destroyed. This is the case also for the neg- 
ative Vsi and there will be ionization levels in the band 
gap also when the vacancy traps a positron. 
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is localized at the vacancy (symmetric Td). 



3. Ionization levels 



The positions of the thermodynamic ionization levels 
for vacancies in semiconductors are also given in Ta- 
ble [H Figure H] shows as an example the determination 
of the level e(-2/-3) for Vca in GaAs without and with 
a trapped positron. The position of the ionization level 
is given by the point where the total energies cross. We 
see that the introduction of the positron lowers the ther- 
modynamical ionization level. According to Table U the 
lowering of the ionization level due to positron trapping is 
a general trend which reflects the lowering of the positron 
energy eigenvalue at the vacancy due to the excess neg- 
ative charge. Slightly surprisingly, the magnitude of the 
lowering is only of the order of 0.1 eV in spite of the 
rather large ion relaxations due to the trapped positron. 
The small change is due to the fact that the magnitudes 
of the ion relaxation in the adjacent charge states before 
or after the positron trapping are rather similar and be- 
cause the excess electron density is rather delocalized and 
does not strongly lower the positron energy eigenvalue. 

The behavior of the thermodynamic ionization levels 
in relation to the positron trapping is an important result 
justifying positron experiments devoted for determining 
ionization levels. Namely, we could think that the neutral 
and negatively charged states could become thermody- 
namically unstable with respect to loosing a bound elec- 
tron when a positron is trapped. However, the lowering 
of the ionization levels indicates that this does not occur. 
Actually, within our model the trapped positron can then 
affect (lower) the measured ionization levels only in the 
case of the level e(0/— ) and maybe also in the case of 
the level e(-/-2). Namely, we expect that only the neu- 
tral or the singly negative charge state (not counting the 
charge of the localized positron) can trap an electron due 
to the positron-induced changes within the positron life- 
time. The more negative charge states effectively repel 
free electrons hindering their trapping. The situation is 
similar to the positron trapping at positively charged va- 
cancieSfS' On the other hand, the ionization level e{+/0) 
between the positive and the neutral charge state is de- 
termined in the positron experiments by the fact that a 



0.7 0.75 0.8 0.85 

(eV) 

FIG. 4: Total energy (zero level arbitrary) of a system com- 
prised of a positron and the Ga vacancy in GaAs as function 
of the charge state of the vacancy and electron chemical po- 
tential fie- The positron is either trapped at the vacancy or 
in the delocalized bulk state. 



positive defect does not trap a positron, and therefore the 
possible change of the vacancy charge state from the neu- 
tral to the negative one does not affect the determination 
of this level. 



4- Positron annihilation characteristics 

Positron lifetimes calculated for the different de- 
fects are also given in Table ID The LDA for the 
electron-positron enhancement effects calculated with 
the Boronski-Nieminen interpolation form underesti- 
mates the positron lifetimes in comparison with exper- 
iment. The underestimation is especially strong for ma- 
terials containing d electrons such as transition metals 
and the III-V compound semiconductors such as GaN. 
Also the lattice constants calculated within the LDA 
for the electron exchange and correlation effects are too 
small compared with the measured ones decreasing the 
positron lifetimes. Therefore rather than the absolute 
lifetime values the differences or the ratios between the 
positron defect and bulk lifetimes are the most important 
figures. The ratios between the defect and bulk lifetimes 
are according to Table U about 1.1-1.25 and 1.5-1.8 for 
typical semiconductors and metals, respectively. For Vca 
in GaN and for Vug the ratios are 1.65 and 1.32, respec- 
tively. The ratios reflect the degree of the localization 
of the positron at the vacancies and, in general, their 
trends are similar to the trends in the positron trapping 
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FIG. 5: The experimental and calculated Doppler spectra 
(ratio to bulk) of triply negative Ga vacancies in GaN (o, 
dash line) and GaAs (•, solid line) . The experimental data 
are from Refs. 153 and respectively. 



energies. It is interesting to note that the change of the 
charge state to a more negative one slightly decreases the 
positron lifetime in typical semiconductors whereas the 
lifetime in Vca in GaN is insensitive to the charge state. 
The deep-level electron vifave functions in Vca in GaN 
are rather delocalized and adding more electrons on the 
deep levels does not appreciably change the total electron 
density and correspondingly the ionic relaxations. 

Figure [5] shows ratio curves between the coincidence 
Doppler broadening momentum distributions of annihi- 
lating electron-positron pairs for vacancies and for the 
corresponding bulk lattices. The measured and calcu- 
lated curves for triply-negative Ga vacancies in GaN and 
in GaAs are shown. The data correspond to the [0001] 
and [001] directions. The calculated distributions are 
convoluted with Gaussian functions with the full width 
at the half maximum of 5.3 and 5.5 xlO^'^ moc for GaN 
and GaAs, respectively, corresponding to the experimen- 
tal resolutions. The calculated curves quantitatively re- 
produce the experimental trends. At low momenta the 
ratio for GaN is higher than that for GaAs reflecting the 
larger reduction of the electron density at the vacancy in 
GaN. At high momenta the GaAs curve is above the GaN 
curve due the contribution of As 3d electrons in GaAs. 
The agreement at high momenta shows that our scheme 
is able to predict the overlap of the positron and core 
electron densities or at least the relative change in the 
positron-core electron overlap between the localized and 
delocalized positron states. 

To enable further studies of trends between different 
materials Table U shows the relative W parameters of 
vacancies, W^d = W^dofcct/Wbuik, that is an experimen- 
tal parameter reflecting the decrease in core annihilation 
compared to the delocalized bulk state. Here Wdcfcct is 
the W parameter (integral over high-momentum part of 
the Doppler spectrum) corresponding to the localized de- 
fect state and Wdofect that of the delocalized bulk state. 
The momentum window used is in all cases 15 — 30 x 
10"'^ moc. The W^d parameter reflects the positron local- 



ization and the extent of the high-momentum core elec- 
trons inside the vacancy. For example, it is interesting 
to note that the VFrei parameter is clearly smaller for Mg 
than for Al. 



IV. DISCUSSION AND COMPARISON WITH 
EXPERIMENT 

Our model gives for vacancies in typical metals and for 
cation vacancies in compound semiconductors positron 
binding energies which are large enough that positron de- 
trapping at elevated temperatures even near the melting 
point is unlikely. The predicted increases in the positron 
lifetimes in trapping are in agreement with experimen- 
tal values. For example, for Al the calculated increase 
is 83 ps in agreement with the experimental one of 85 
ps (Ref. 112) and for for Vca in GaN the calculated and 
measured increases are 85 and 75 ps (Ref. HO), respec- 
tively. The calculated and measured electron-positron 
momentum distributions at high momenta also show a 
good correspondence (Fig. Ogives an example) meaning 
that our model is able to reproduce the positron overlap 
with ion cores in a reasonable manner. 

The agreement between theory and experiment be- 
comes less clear when the predicted positron trapping 
energy decreases. In the case of positron trapping at va- 
cancies in Mg the calculated positron lifetime increase is 
from 219 ps to 289 ps whereas the much smaller increase 
from 225 ps to 255 ps has been measured between well a 
annealed sample and a sample with thermally generated 
vacancies?^ The calculated large lifetime increase is also 
reflected in a rather small Wrci parameter of 0.56 for ^Mg- 
However, the calculated positron trapping energy of 0.34 
eV is in agreement with the estimate of 0.3 — 0.4 eV by 
Hautojarvi et al^ One possible source of the theory- 
experiment disagreement could be difficulties in extract- 
ing bulk and vacancy lifetimes for Mg in measurements. 

Our calculations predict a vanishingly small positron 
trapping energy at vacancies in Si, of the order of ther- 
mal energy at room temperature. The trapped positron 
state does not exist when the vacancy is relaxed without 
the influence of the localized positron, which in princi- 
ple prevents positron trapping at the vacancy at zero 
temperature. Our result is, however, even qualitatively 
wrong since the predicted trapping energy is negative. 
However, the behavior of Vsi during positron trapping is 
unique reflecting the very flat energy landscape the ions 
feel around the vacancy. The flatness is due to the dif- 
ferent competing possibilities for bonding and rebonding 
in a covalent material. The flat energy landscape is also 
behind the scatter of the DFT results for the structure 
and energetics of the Si vacancy. Only recently, when 
calculations with very large supercells (up to 1000 Si 
atoms) have become possible the results show a satis- 
factory numerical convergence!^ Figure [Ijb) shows that 
the energy landscape of the Si vacancy with a trapped 
positron is even flatter than that of the plain vacancy. 
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This means, as discussed above, that the entropy con- 
tribution should be taken into account when describing 
the trapped state at finite temperatures. However, one 
should bear in mind that the errors arising, e.g., from 
the LDA's for the electron-positron correlation energy 
and for the electron-electron exchange and correlation 
may be of the order of tenths of an eV. For example, 
the scatter in the calculated formation energies for Vsi is 
of this order or even larger. Thus, our scheme may de- 
scribe even qualitatively incorrectly the actual positron 
trapping process in borderline cases such as Vsi. 

The flat energy landscape for Vsi would introduce 
strong temperature-dependence to positron trapping and 
dctrapping processes. In our model [Fig. [Ijb)] the 
positron trapping would, in principle, be possible at finite 
temperatures because part of the time the vacancy vol- 
ume is so large that a bound positron state exists. This 
probability is, however, strongly temperature-dependent 
and in the experiments the trapping rate as a function 
of temperature is seen to be rather constant except for 
the case of negative defects for which the trapping rate 
decreases with increasing temperaturej^ 

The concept of positron binding energy seems to be a 
difficult one to define or at least it is difficult to get a 
quantitative agreement between experiment and theory 
even if the computational results were exact. First of 
all, one of the assumptions behind the relation between 
trapping and dctrapping rates^^ typically used in the in- 
terpretation of experimental data is that the excitations 
of the positron are decoupled from the excitations of the 
system. According to our calculation this clearly is not 
the case. Secondly, the positron dctrapping energy (the 
threshold energy needed to detrap a localized positron) 
in the case of Vsi strongly depends on the dctrapping 
mechanism (fast vertical transition vs slow transition due 
to phonons) and on the ionic structure of the defect at 
the instant of dctrapping [see Fig. [T]^b)]. In the case of 
this kind of an energy landscape it is difficult to draw 
conclusions about the actual positron trapping and de- 
trapping processes. According to our calculations the 
strain energy stored in the relaxation of the vacancy is 
so large that all the released energy can be stored in it 
(analogously to the trapping of a small polaron). How- 
ever, although the energy is conserved between the ini- 
tial and final states the differing time scales in electronic 
(transition of the positron from delocalized to localized 
state) and ionic processes (relaxation of the vacancy dur- 
ing positron trapping) complicate the picture. Although 
our results may not in all cases be even qualitatively cor- 
rect they clearly suggest that the models for positron 
trapping in semiconductors- need refinement. 

Recent measurements of heavily As-doped Si indicate 
that positron trapping energies at vacancy defects in 
highly As-doped Si may be low and thermal dctrapping is 
possible. Namely Kuitunen et al. found that positron 
dctrapping happens from the Si vacancy decorated by 
three As atoms (Vsi — Asa) at temperatures above 500 K. 
For vacancies decorated with one or two As atoms they 



did not notice positron dctrapping. Using the relation 
between the dctrapping and trapping rates derived by 
Manninen and Nieminen^ Kuitunen et al. determined 
for Vsi — As3 the trapping energy of 0.27 eV. Our scheme 
gives a clearly smaller positron trapping energy of 0.06- 
0.1 eV for Visi — Asa. Similarly to the case of Vsi we do not 
find for Vsi — As and Vsi — As2 an energetically favored 
trapped positron state but just a metastable configura- 
tion in disagreement with the experimental trend. The 
trend in our trapping energies is exactly the opposite; the 
trapping energy increases with the increasing number of 
As atoms (n = 0, 3). This is because the large As ions 
around the vacancy do not relax inwards as strongly as 
the neighboring Si ions so that the energy stored in the 
ionic lattice in the positron trapping process is smaller for 
the As decorated vacancy than for the clean vacancy. A 
well-known shallow positron trap in Si with open volume 
is the complex formed by a vacancy and an interstitial 
oxygen. For it, the reported positron binding energies 
are of the order of 40-50 meV (Refs. Ulii). The open 
volume of the V — complex is so small (see Ref . \at\ \ that 
we do not expect it to trap positrons in our calculations. 

The fact that a bound positron state can be found at 
Vgc or at Vas in GaAs when freezing the ions at ideal 
positions means that the bound positron states are very 
close to appear, and an improvement in the theoretical 
description could lead to bound positron states also for 
optimized ion positions. In experimental works^2ii2ii^ 
positron lifetime components between 279 and 292 ps 
are assigned to Vce- The measured lifetimes can be con- 
trasted to the measured bulk Hfetime of 228 ps (Ref. lish. 
These vacancy lifetimes are already quite close to the the- 
oretical estimate of 316 ps for an ideal divacancy in Ge 
(the corresponding bulk lifetime is 229 ps)pi2. For an ideal 
neutral monovacancy we get the lifetime of 246 ps which 
is only 33 ps longer than our bulk lifetime. These com- 
parisons suggest that the experimentally observed life- 
time components are too long to be explained by annihi- 
lation at ideal monovacancy-size defects. Measurementa^ 
for n-type GaAs show the lifetimes of 257 ps and 295 ps 
(the experimental bulk lifetime is 231 ps). These were 
assigned to negative and neutral As vacancies possibly 
associated with impurity atoms on the basis of measured 
ionization levels and corresponding old computational re- 
sults^i^ which suggested that the levels of Vas are near 
the conduction band. The measurements gave ionization 
levels of e(+/0) = ~ 0.140 eV and e(0/-) = - 
0.030 eV, where Ec is the conductance band minimum. 
Our first-principles results, which are in accordance with 
those in Ref. [53, indicate that these ionization levels of 
clean Vas in GaAs are close to the top of the valence 
band. This supports also the conclusions that the de- 
fects observed in positron measurements are not clean 
vacancies. 

Finally, our calculations for in GaN suggest that 
bound positron states at small anion vacancies in com- 
pound semiconductors are not possible. This should be 
contrasted with a positron annihilation study in which 
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the detection of N vacancies is reported.-^ In the inter- 
pretation, however, the short Ufetime component is asso- 
ciated with VN-impurity complexes and not with isolated 

V. CONCLUSIONS 

We have studied, using first-principles calculations, the 
role of lattice relaxations around vacancies in the positron 
trapping process in various metals and semiconductors. 
In metals the trapping energies are typically of the or- 
der of 1. . . 2 eV. The lattice relaxes due to the trapped 
positron and the positron annihilation parameters change 
but especially the qualitative effects are small. The most 
important difference between typical metals and semi- 
conductors is the magnitude of the strain energy stored 
in the lattice compared to the lowering of the positron 
energy eigenvalue. For semiconductors these two can be 
of the same order of magnitude because of the smaller 
reduction of nucleus repulsion in the trapping process. 
This leads to rather small values for the trapping energy. 

In the case of covalently bonded elemental semicon- 
ductors such as Si and Ge the energy landscape of the 
positron- vacancy system is extremely flat which suggests 
that entropic considerations have to be taken into ac- 
count when discussing the ionic structure of the vacancy 
with a localized positron at finite temperatures. For Si 
the calculations give, in disagreement with experiments, 
vanishingly small or rather sightly negative trapping en- 
ergies suggesting that trapping of thermalized positrons 
would not even be energetically favorable. 

Also the vacancy in Ge and anion vacancies in com- 
pound semiconductors are challenging tests for theoret- 
ical methods since in the calculations there is either no 



bound positron state at the vacancy or the trapping is en- 
ergetically unfavorable and there is no local energy min- 
imum configuration at which the positron is trapped at 
the vacancy. Furthermore, for the vacancy in Si the con- 
figuration with the trapped positron is just a metastable 
state while in the ground state the positron is in the 
delocalized bulk state. An interesting finding is that a 
localized positron cancels all the Jahn- Teller distortions 
we observed when having no positron at the vacancy. 

In general, the higher the predicted positron trapping 
energy is the better is the agreement between our results 
and experiments. However, also in the case of defects in 
Si we get a surprisingly good agreement in calculated life- 
times and Doppler broadening spectra for the metastable 
state. 

We have also studied the effect of the localized positron 
on the electronic structure of the vacancies in semicon- 
ductors by evaluating thermodynamical ionization levels 
of vacancies. The positron-induced changes are usually 
only of the order of 0.1 eV. In general, the levels move 
closer to the valence band maximum. A positron-induced 
change in the defect charge state in this model is likely 
only in the case of a neutral defect since positive ones do 
not trap positrons and negative ones electrons within the 
lifetime of the trapped positron. 
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